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Abstract 

The expected large mass rrit of the top quark provides for the first time a chance to discuss 
the bound-states of the corresponding quantum system as the nonabehan generahzation 
of positronium, using the full power of relativistic quantum field theoretic methods which 
are available for weakly bound systems. Thus our approach differs in principle from the 
one used in the vast phenomenological literature on quarkonium potentials. We emphasize 
especially the corrections of energy levels which are of order a'^rrit or numerically com- 
parable to that order, and which have no counterpart in the 'relativistic' corrections of 
QED. In contrast to previous computations we give analytic expressions for all contribu- 
tions considered in our present work, hopefully preparing the ground for further similar 
calculations. 
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1 Introduction 



Perturbative expansions in the coupling constant in quantum field theory possess two 
types of applications, the calculation of scattering processes and the computation of pro- 
cesses involving weakly bound systems. Many of the successes of quantum electrodynamics 
(QED) are, in fact, related to positronium, i.e. to the second one of the aforementioned 
applications. The proper starting point for any bound-state calculation in quantum field 
theory is an integral equation, comprising an infinite sum of Feynman graphs. The Bethe- 
Salpeter (BS) equation 0] fulfills this task and it is well known that in the limit of binding 
energies of O(a^m) the Schrodinger equation with static Coulomb attraction is obtained. 
The computation of higher order corrections to the Bohr-levels, however, turned out to 
be far from trivial. It was recognized, though, relatively late that, at least conceptually, 
substantial progress with respect to a systematic treatment results from a consistent use 
of a perturbation theory geared to the original BS equation 0. In that manner, at the 
same time, nonrelativistic expansions as implied by Hamiltonian approaches with suc- 
cessive Fouldy-Wouthuysen transformations are avoided. Within the BS-technique, 
however, it is desirable to have an exactly solvable zero order equation different from the 
Schrodinger equation, otherwise e.g. the approximation procedure for the wave function 
lacks sufficient transparency, especially in higher orders. One of the advantages of the BS 
approach to perturbation theory is the freedom to select a different zero order equation. 
Of course, in that case, certain corrections included already at the zero level are to be 
properly subtracted out in higher orders. An especially useful zero order equation has 
been proposed some time ago by Barbieri and Remiddi (BR equation [^). Still, one of 
the most annoying features of all bound-state calculations remains the pivotal rule played 
by the Coulomb gauge. In other gauges, e.g. already the (in QED vanishing) corrections 
0{a^m) of the Bohr levels imply to take into account an infinte set of Feynman graphs 
Only in very special cases, when certain subsets of graphs can be shown to represent 
together a gauge-independent correction, another more suitable gauge may be choosen. 
By contrast to QED the vast literature on bound state problems in quantum chromody- 
namics (QCD) adheres to a description of the quark- antiquark system by the Schrodinger 
equation with corrections 'motivated' by QCD |^. As long as a relatively small number 
of parameters suffices for an adequate phenomenological description of observed quantum 
levels, this approach undoubtably has an ample practical justification. However, again and 
again certain deviations from such phenomenological description are reported 0. Thus 
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also for this reason a return to more rigorous QCD arguments remains as desirable as 
ever. The standard literature on this subject almost exclusively is based on nonrelativis- 
tic expansions or on the calculation of purely static Coulomb forces |^ . Also very often 
potentials with higher order corrections as determined from on-shell quarkonia scattering 
are used 0. In that cases relevant off-shell effects may be even lost which are typical for 
higher order corrections. On the other hand, from the point of view of relativistic quantum 
field theory as elaborated in the abelian case of positronium, a similar, more systematic 
approach seems desirable, the more so because the basic techniques are well developed. In 
addition, at least in one case, namely the decay of S-wave quarkonium, the result of a full 
BS-perturbation calculation |jlO[ , including the QCD corrections to the bound state wave 
function, yields a result very different from the one which considered only the corrections 



to the quark antiquark annihilation alone ITI . 

In this connection the relatively large size of the running coupling constant even at 
high energies represents a well known problem, together with large coefficients from a 
perturbative expansion. For this reason e.g. problems arise in the comparison of the cou- 
pling constant as determined from scattering experiments within the minimal subtraction 



scheme (MS), with the coupling constant to be used in a consistent weak bound-state 
approach. The philosophy within our present work will be that the orders of magnitude, as 
determined from ajjg will be used for estimates, but that we shall imply a determination 
of as by some physical observable (e.g. energy levels, cf. the remarks after eq. (|4.31| ) ) of 
the quarkonium system itself. In that way delicate correlations of 'genuine' orders of 
from basically different types of experiments are avoided. 

Of course, the quarkonium system also differs profoundly from positronium because 
of the confinement of quarks and gluons. However, the phenomenological success of the 
nonrelativistic quarkonium model can be explained by the fact that the bound states of 
heavy quarkonia are deep in the Coulomb funnel and thus sufficiently far away from the 
confinement part of the potential. Estimates in the early 80-s [|T^ of nonperturbative ef- 
fects, describing the 'tail' of confinement by a gluon condensate |T3[ suggested that only 
with quark masses well above about 50GeV the importance of such nonperturbative ef- 
fects for low Bohr quantum numbers may decrease sufficiently to make perturbative 'field 
theoretical' level corrections competitive and observable. From high precision electro- weak 
experiments of the LEP collaborations, the mass range of the top quark now seems to 
be etablished to lie in the range lOO-lSOGey WM. Thus for the first time a nonabelian 
bound state quarkonium system seems to fullfil the high mass requirement for a genuine 
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field theoretic approach. Unfortunately the drawback of this situatuion is that the weak 
decay t ^ b + W broadens the energy levels |15| for increasing mass rrit so that above 
nit ~ 14:0GeV individual levels effectively disappear. Even for nit ~ lOOGe^ this broad- 
ening of energy levels certainly already makes the resolution of different Bohr levels of 
0{ma'^ /An"^) ^ IGeV (for a ^ 0.2) difficult to observe. The first field theoretic correc- 
tion for nonabelian QCD already starts at O(ma^) 16|, and the corresponding shift 



depends on the principal and the angular momentum quantum number. Although one 
thus has to fully acknowledge the problems related to experimental observations of cor- 
rections to O(a^m) for the energy levels of the toponium system, we have been tempted 
to trust the ingenuity of experimentalists to eventually arrive at sufficiently precise data 



I7| , together with the hope that mt does not lie too far above lOOGeV^. In that case, for 



the first time, genuine field-theoretic consequences of QCD could be tested at something 
like a nonabelian positronium system. Of course, the well known 'relativistic' correc- 
tions, corresponding to the same type of graphs as in the abelian case are present here 
as well. However, already considering only vertex corrections, a difference to the abelian 
case at O(a^m) was discovered a long time ago |16|. In addition typical other nonabelian 
contributions may appear at this level as well. In contrast to scattering processes, the de- 
termination of the order to which a certain graph contributes in relativistic bound state 
perturbation theory in each case requires a special analysis. Other terms (from QED 
and weak interaction), incidentally, may be of the same numerical order O(a^m). The 
Coulomb gauge also entails peculiar additional nonlocal interaction terms in the effective 



action, appearing e.g. in the path integral formulation |T8[. We also investigate the effect 
of those terms here. 

In section 2 we recall some basic facts about BS-perturbation theory and about the 
BR equation for nonabelian weakly bound onium-systems. 

Tree graphs leading among others to the well known 0(Q;^m) corrections are discussed 
in sect. 3. As indicated already above, the one loop vacuum polarization (sect. 4.1) provides 
a correction term O(a^m) in the nonabelian case. Here we take the opportunity to point 
out that effects from some of the lighter quarks (charm, bottom) in the toponium system 
must be treated more carefully than it is usually done by including the quarks only in the 
number of flavours of a running coupling constant. 

In sect. 4.2 we revisit the old calculation of Duncan [|l6l within our present formalism, 
avoiding some approximations made in this early computation which allows us to make 
even some statements on the term 0{a^m). Finally sect. 5 is devoted to an exploration of 
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possible further corrections of the same numerical order of magnitude as 0{a'^m). We list 
several candidates of relevant QCD graphs and we show in cases of most simple two loop 
graphs that such corrections may well be relevant. Beside these graphs from QCD, also 
corrections from the weak interaction and QED may turn out to contribute to this order, 
but the nonlocal Schwinger-Christ-Lee type graphs, peculiar to the Coulomb gauge, are 
irrelevant to 0(a;^m) as shown by an explicit calculation. 



2 BS-Perturbation of the BR Equation 

A correct formulation of QCD in Coulomb gauge entails not only Faddeev-Popov-ghost 
terms but also the inclusion of nonlocal interaction terms |T^. Therefore, the full La- 
grangian reads (a=l,...,8 for SU(3)): 

£ = -lF;^F'^'^^' + j2^^i^^D-m,)^J + B'^{^,A])-r]'^^\5a0^+9fab^^^^ (2.1) 



where the Lagrange multiplier B"^ guarantees the Coulomb gauge, and where 

= d.-tgT'^A;, (2.2) 
F;, = d,At-d,A';^ + gUcAlA:. (2.3) 



f 1 and f 2 are given in [|T8[ and are discussed more explicitely below. The above Lagrangian 
will include all effects of the strong interaction, but, as we will show, QED and weak 
corrections may also give a contribution within the numerical order of our main interest 
(0(a^)). 

In our notation the BS equation in terms of Feynman amplitudes for K and S reads 

as 

Xf (p;P) = ~^S^^iJ+p)S,,,i-^+p) I i^^ir,,^.,,„^.„(P,p,p')x^'?"(p';^'). (2.4) 

X denotes the BS wave function, S the exact fermion propagators, and K is the sum of 
all two fermion irreducible graphs. Furthermore, we have introduced relative momenta p 
and p', a total momentum P = pi — p2, and we choose a frame where P = {Pq, 0). 



The notation can be read off the pictorial representation in fig.l. i,j are collective 
indices for spin (noted a, p) and colour (noted a,P). It is well known that the dominant 
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part in K for weak binding (a — 0) is the one-Coulomb gluon exchange which results 
in an ordinary Schrodinger equation with static Coulomb potential. This result is even 
independent of the chosen gauge in the ladder approximation by a simple scaling argument 
Pq ^ Oima^), \p\ ^ 0{ma), Pq ^ 2m — 0{ma^) p9|] : 

K^K^:= - (-^7^7,% (2.5) 

In equation ( f^.Sj ) we have already used the fact that only colour singlet states can form 
bound states because the Coulomb potential is repulsive for colour octets. The colour 
trace will always be understood to be already done, leading to the definition 

a = = -as- (2.01 

3 47r 3 ^ ^ 

in terms of the usual strong coupling constant as- Because the above mentioned nonrela- 
tivistic limit of the BR equation contains the projection operators A^, defined below, it is 
awkward to calculate the so-called relativistic corrections in a straightforward way within 
the framework of BS perturbation theory, starting from ( ^.41 ) with ( |2.5| ). Therefore, we 
use the BR equation instead of the Schrodinger equation. It is obtained from ( |2.4| ) by 
the substitutions 

S 



1 



7P — m 

K ^ Kbr=[ioA+\^AX',A^'^'^'1o]p',p"K, (2.7) 



with the projection operators 



Ep + dp ± Pm 



= ^(l±7o) 



k . -^^^ 2E, 2 2 _ 



and 

K ^ - 

{p-p'Y Ep + mEp> + m ^ + 2Ep ^ Pq + 2E. 

m denotes the mass of the (heavy) quark. The solutions of the equations obtained in this 
manner are exact. The (colour singlet, normalized [Q) eigenfunctions read {uj := Ep — ^) 
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x{p) 



x{p) 



A+A+rA-A- Ej2Ep + Po{Po - 2E, 



/I^{Ep + m) 



{po - LJ + ie){po + uj - ie) 



7o(A+A+rA-A-)*7o EpJ2Ep + Po{Po - 2Ep 



-(j){p) 



(2.10) 



{po-uj + ie) {po + uj- ie) 



^{Ep + m) 



and belong to the spectrum of bound states 



Po = MO = 2m^l- — ^2m-m—^ + 0(a«). (2.11) 

In eqs. ( p.lO| ) F is a constant 4x4 matrix which represents the spin state of the particle- 
antiparticle system: 

r = {?^' ■ (2.12) 

is simply the normalized solution of the Schrodinger equation in momentum space, 
depending on the usual quantum numbers (n, /, m) |]2T|, a±i, oq describe the triplet states. 
In the following it will often be sufficient to use the nonrelativistic approximations of eqs. 

X(P)"^ = , "^T^ . 0(p)r, (2.13) 
Po — 1^ + te 

X(P)"^ = , . <l>*m-loT*lo). (2.14) 

Pq — uj"^ + te 

Perturbation theory for the BS equation starts from an exactly solvable equation, in 
our case the BR equation for the Green function Gq of the scattering of two fermions 

tGo = -Do + DoKoGo. (2.15) 

Dq is the product of two zero order propagators, Kq the corresponding kernel. The exact 
Green function may be represented as 

1 oo 

G = i: Xffp 7;r^rf + Ore, = Go E(^G'o)^ (2.16) 
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where the corrections are contained in the insertions H. Bound state poles M„ contribute, 
of course, only for Pq < 2m. It is easy to show how H can be expressed by the full kernel 
K and the full propagators D: 

H = -K + Ko + iD-^ -iDq\ (2.17) 

Since the corrections to the external propagators contribute only to O(a^) |2^, the per- 
turbation kernel is essentially the negative difference of the exact BS-kernel and of the 
zero order approximation. 

Expanding both sides of equation (|2.16|) in powers of Pq ~ with M° from (|2.11|) 



the mass shift is obtained 0: 

AM = (/io)(l + (/ii)) + {hogih^) + Oih"). (2.18) 
Here the BS-expectation values are defined as e.g. 

We emphasize the four-dimensional p-integrations which correspond to the generic case, 
rather than the usual three dimensional ones in a completely nonrelativistic expansion. 
In ( |2.18|) hi and Qi represent the expansion coefficients of H and Gq, respectively, i.e. 



H = J2hn{Po-My (2.20) 

n=0 

oo 

Go = Y.9niPo-My'' (2.21) 

n=0 

3 QCD Tree Diagrams 

The contributions stemming from the tree diagrams 2. a to 2.c are well-known from the 
abelian case. Fig. 2. a is peculiar for the use of a different zero order equation than the 
Schrodinger equation. It contains the difference between the exact one Coulomb-gluon 
exchange and the BR-Kernel ( p.7|) . The exchange of one transverse gluon is represented 
by graph 2.b, and fig. 2.c shows the annihilation graph. The latter does not contribute in 
our nonabelian case. 

For the sake of completeness and in order to illustrate the present formalism, we exhibit 
first the results for the tree graphs as well. The perturbation kernel for the Coulomb-gluon 
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exchange 



(3.22) 



is needed for the calculation of the energy shift induced by fig. 2. a using Eqs. ( ^.181 ), 
( CT ), (EUD and to For the spin-singlet we have: 



AMi., 



He + Kbr) 

2 



a 



ma 



Tia f ay f 



pp 



5, 



10 



+ 



Pom J {2Tif J (27r)3' 
1 



{p — p'y 



1 



+ 0(a^). 



16n4 16n3(/ + l/2) 
The contribution from the transverse gluon (fig. l.b) 

■TT -A i k ^ / . 'ij'ik- 



(3.23) 



(3.24) 



with 



q = p - p, 



(3.25) 



gives rise to a spin singlet-triplet (magnetic hyperfine) splitting. Because of the 7-^ matri- 
ces, the projectors from both wave functions annihilate (|3.24|) . This means that two 
factors fry, contained in A^, are needed for a nonzero result which in turn gives rise to 
two extra orders of a. By this mechanism we arrive at the well known contribution 0{a^) 
from this graph. For the spin-singlet the mass shift reads 



AMi. 



2TTa , 



b,S=0 



1^ 



m^ 



+ 2 



4/ 

ma [ 



^10 



d^p 
(2^ 
3 



<'V^.,-,,,Mp)_§)^(p-)]+0(a») 



8n4 8n3 16n3(/ + | 



(27r)3 

+ Oia^). 



(3.26) 



The evaluation of the singlet-triplet splitting requires some awkward Dirac-algebra, but 
the final result may be brought in a quite transparent form (where one recognizes this 
expression as the well known spin-spin and spin-orbit interaction, adapted to the present 
problem, cf. e.g. ||2| 



AM 



ortho,m 



AM, 



para 



2'Ka 

m^ 
2Txa A 
~m?^^ 



d^p / rfV0*(p')(l + 



5(F) + 



1 f 2 


- 3 ra^ P 


3i{fxp){a*^ X 


An 


1 T^l ^ 


Attt 3 



(3.27) 
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For our purpose it is sufficient to calculate the singlet-triplet splitting for the ground state 
given by the term proportional to 6{r): 

AM,,,,„ - AM,ara = ^l*(0)r = (3-28) 

It should be mentioned that this result differs from the positronium case, because the 
graph 2.C gives no contribution here. Moreover, we want to emphasize again that according 
to our 'purist' field theoretic point of view, integrals like ( p.27| ) and (|3.28|) , of course, are 
not to be evaluated between phenomenological wave-functions. 



4 One Loop Corrections 

4.1 One Loop Vacuum Polarization 

In the case of positronium no massless particles can contribute to vacuum polarization and 
so this effect is only of order of magnitude a^. In contrast to this, QCD contains massless 
gluons and light quark flavours which may contribute significantly to the spectrum. 

Using standard techniques for evaluating the corresponding integrals one obtains for 
the loop in fig. 3.a (A^=3 for QCD): 



(27r)^r2(g-f)2 



'^^'^ f\p^\^f + L-2\n2- (4.29) 



3^2 - L- .3 

7 c 

-e( 2 In 2 - 7 In Ati) In g ^ + - In^ g ^ + e ■ const + 0(e^)] 

3 2 



with 



1 A- D 

V = --7 + ln47r, e = — ^ (4.30) 

We have written down this result correct to orders €■ f{q) in order to make it applicable in 
the two loop calculation. The graph 3.b is more difficult to calculate, because it contains 
go terms. These terms can be completely avoided if we carry out the pq integrations ffist 
(cf. ( p^.l9D and ( [^.18D ). The result can be expanded in powers of p and a in order to show 



that the effect of go is of 0{a^). With this simplification graph 3.b is exactly calculable: 
^tb = ^^^m'{V - \nf) + 16(7 - 81n2)g-'2] + 0(e) (4.31) 
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Our renormalization prescription consists of a subtraction at the point q = (0, fl), where /i 
has to be of the order am to avoid large logarithmic contributions from higher orders. This 
seems to be the natural renormalization prescription for bound state problems, because 
also in the BS expectation values ( |2.19| ), the Bohr momentum am together with ^ 
0{a'^m) provides the dominant parts of the integrals. 

After renormalization, the contribution from the gluonic vacuum polarization (with 



the colour trace already done) reads |jT6 



-iHg = -i-fo®loHg (4.32) 

The expectation value of this expression can be obtained by performing the Fourier trans- 
formation into coordinate space, where the integrations can be done analytically (see 
Appendix A). Our surprisingly simple result is 



AM, = (Hg) = -ma^^^[^,{n + / + 1) + 7 + In i^] + ©(a^) (4.33) 

16?™^ am 



where \E'„ is the n-th logarithmic derivation of the gamma function and 7 denotes Euler's 
constant. The closed form of Eq. ([4.33|) was not obtained in previous calculations. 



Now we turn to the contribution from the fermion loops fig. 3.c. In the literature the 
lighter quarks are usually taken as massless (and 'absorbed' in the number of flavours 



appearing in the running coupling constant) or even ignored [16] [Q, but we will show 
that they do contribute within the order of interest and, furthermore, the explicit depen- 
dence on the masses of the lighter quarks is important. As pointed out already in the 
introduction, this is due to the fact that the top quark is expected to lie above lOOGeV 



2J] and therefore the bottom and charm quark can neither be taken as relatively massless 
nor as relatively super-heavy as compared to the natural mass scale am. 

The finite part of the self energy in fig. 3.c is a well known quantity p5[ for arbitary 
masses mj of the quark: 

4^2 ^ '-J 



^i. = -^g^/(g^m?). (4.34) 



We approximate in the exact solution 

/((f^,m^) = / dxx{l — x)\n[x{l — x)!]"^ + 
■' Jo 



m 



llnm^2-^ + /(p), (4.35) 
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/(p) 



2p 1 



3+-(l-2,)Vl+4pln.^^_^. 



2 



for later convenience /(p) by 



(4.36) 
(4.37) 

(4.38) 



This agrees with the original /(p) better than 1% within the whole integration region. 

It seems instructive to transform into coordinate space in order to obtain the potential, 
effectively produced by this fermionic vacuum polarization: 



a 



5 

[Ei(— rmjcs^ 



5 1 , , p 5, 

6 2 ^m? ' 



(4.39) 



The mass shift can be obtained from ( |4.39| ) in closed form using the integral formula 



r e-^''Ei{-ax)dx = -^[e-^"Ei(-ax) + ln(l + -) - 
Jo p a 



-Ei{-{a + p)x)]. 



(4.40) 



Thus a useful expression for the energy shift induced by fermionic vacuum polarization 
with arbitrary masses rrij reads 



AM], 



ma 



3 2n-2l-2 



k=0 



d 



2l+l+k 



— [- \n{l + (3 



am 



dp 



nmica 



(5^1 



5 1, ,p' 

6 2 mi 



with 



{n-l-l)\ ( d 



■ A;![(n + /)!]3 \dp; 
For states up to n = 3 we write this result more explicitly as 



p=0 



(4.41) 



(4.42) 



ma 



^ nmj . 



+ 63)2 



am 



+ ln- 



66 + 



(4.43) 



with A„j from Tab. 2, 
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using the shorthand 



n 


/ 


^ am ' 


1 







a 




2 







— + a 




2 


1 




+ + a 




3 









+ a 


3 


1 


a' 


- + + 


+ a 


3 


2 


5 


+ + + ia^ 


+ a 



Tab. 2 



-1 ._ 



1 + 



am 



(4.44) 



Only for mj » am this gives an 0{a^) Uehhng term, modified by off-shell subtraction, 
but for mj ^ it becomes an 0{a^) contribution, which means that eq. ( [4.43| ) interpolates 



numerically in a range of two orders in a. Therefore ( ^.43| ) must be definitely taken into 
account for quarks with m^ am at 0{a^). 



5 Vertex Corrections 

The one loop corrections to the vertex together with self-energy insertions into the fermion 
lines (fig. 3.d) in the abelian case (positronium), are known to provide corrections only 
of 0{a^). The reason for this is the "classical" Ward identity which continues to relate 
those contributions in such a way that the sum of these terms vanishes at |g| 0. 
This Ward identity happens to continue to hold even in the Coulomb gauge and even 
in the nonabelian case 0, but only for the vertex corrections referring to the Coulomb 
component of the gauge field. However, the presence of the gluon splitting graphs 3.e and 
3. f produces a contribution already to 0{a^) |jTB[. A simple dispersion theoretic argument 
allows to understand this difference: In the abelian case the first graph 3.d in the variable 
|g| for go = has a cut for Re|g| < 2m. Thus corrections in g, e.g. in the electron 
form factor Fi, for symmetry reasons must be of order g^, because Fi(g^) is regular at 
g ^ 0. This is no longer the case with the mass zero intermediate state allowed in 3.e 
and 3.f. The first - and to our knowledge only- computation of the nonabelian vertex 
corrections in the sense of our present approach was performed in ref. This work 
contains certain approximations which we wanted to avoid in order to prepare the ground 
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for future calculations even at the level 0{a^). We thus make a systematic expansion 
and solve the remaining integrals analytically which contain contributions of the order of 
interest. The vertex correction of fig. 3.e after performing the colour trace becomes 

x[/^ 7° -r^^ 7l 

7P + 7r — m 7p — ^yr — m 

= ^^^^(7V.(1,J>)7°-7V.(-1,P')71 (5.1) 
q ^ 



where ( e = ±1) 



Vki{e,p) := TTT^^J^ — ^-4* + — 7 — ^ ^ 5.2 



After the tq integration it proves useful to proceed with the po integrations (cf. Eq. ( p. 19 



) ), where in contrast to ref. |T6[, who approximates already at this point, we took into 
account also the pole arising from the denominator of Eq. ( ^.2|) . This results in 

vUe.p) = -^/ (2t)3~(F-^g-)T -^b>1 (5-3) 



(5.4) 



The calculation of the 7 trace to lowest order requires the inclusion of the 7p-terms in 
the wave functions ( |2.10|) stemming from the projection operators ( |2.8| ). After performing 
the 7 trace we expand in terms of p and ma^ which enables us to combine both terms in 
eq. (O): 

367rV r d^r - ^ 1 



= J (2^)3 (f-qy 2r-E, + ^^^^^^ ^'-'^ 

Now Eq. ( p. 51 ) may be evaluated exactly in terms of dilogarithms and the result has a cut 
for |g| < 0, but no pole at |g| = 0. It can be formally expanded for |g| > to 0(|g|): 

H3. = |^(^-f + 0(,-)). (5.6) 

\q\m 8 3 



The first term in this expansion has been obtained in |TB|, the second one is the expected 
contribution to 0{a^). The BS - expectation value of ( p.6| ) becomes to 0{a'^): 

^ , ^ Qma^ 1 
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It remains to calculate the vertex corrections with two transverse gluons , graph 3.f. 
At first sight it seems that this graph would give a contribution to order because no 
fry terms are needed from the wave functions. This, as we will see, is not true because 
the leading (constant) term vanishes as a consequence of renormalization and accordingly 
graph 3.f has been estimated to be of order O(a^) in ref. |]l6l. Here we use an approach 
which explicitely provides at least part of the exact result of this contribution. The vertex 
part of the graph 3.f reads: 



^ _ jqi - ri){qk -Tk) _ Tifk ^ {qi - ri)r{q- f)rfc \ 
y [q — r)"^ r^{q — r)"^ j 

With the gamma-trace to relative 0{a'^) and using Feynman parametrization, the effective 
vertex part becomes 

- //^ td. r ci.(L-i-3.WL-|)(l+ 

(5.9) 

where 

L = (yqo — xpoY + k"^ + 2xpk — 2yqk — x{p^ — m^) — yq'^. (5.10) 

For our estimate it is sufficient to consider the part stemming from the constant term in 
the second bracket in Eq. ( |5.9| ). The k integration can be done easily leaving us with a 
finite part 

where 

M = x^m^-w, (5.12) 
w = y{l — y)(f + x{l — x){p^ — m^) + 2xypq. (5.13) 

This expression cannot be expanded in terms of w because this would yield a spurious 
linear divergence in the next order from the q^ term which would indicate an equally 
spurious 0(a;^) contribution. Therefore, we expand in terms of {w + y{l ~y)q'^) and solve 
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the leading part analytically in terms of dilogarithms. Expanding the result in terms of 
(f^, one has 

Vlif = ^^{2-4[l + M|^)]^} + otherO(a^) (5.14) 

Since ( p.l4| ) does not contains a term oc the vertex correction due to two transverse 
gluons does not contribute to O(a^). 



5.1 Other QCD graphs 

To one loop order also the box graph 3.g occurs in the correction to the BR kernel. 
It possesses an exact counterpart in QED and is known to contribute only to O(a^) 



27| . Furthermore we have also investigated the two- loop vertex-correction depicted in 



fig. 3.h. Of course this correction is but one of several two loop vertex corrections. The 
renormalization must take into account the whole set of this graphs. Nevertheless, it seems 
that this graph, after proper renormalization, yields an O(a^) contribution. 



6 Other Corrections 

6.1 Two Loop Vacuum Polarization 

As pointed out already in subsection 4.1, the usual renormalization group arguments re- 
lying on massless quarks in the running coupling constant do not consistently include 
the effect of 'realistic' quark masses in the toponium system, when a systematic BS per- 
turbation is attempted. However, in the one loop case finite quark masses gave terms of 
numerical order 0{a^), therefore the same may be expected here, leading to corrections of 
0{a^). On the other hand, in a full calculation of effects of O(a^) two loops with gluons 
cannot be neglected. Altough it is enough to consider the two loop vacuum polarization 
for Coulomb gluons only, the computation of all those graphs is beyond the scope of our 
present paper. We just want to indicate how already the graphs 4.a-4.c yield contributions 
of this order. The result of the self energy part from 4. a to 4.c are given in tab. 3. Let us 
start with graph 4. a. Performing the zero component integrations of momenta results in 
(including a symmetry factor 1/2) 

-iH,,, = -z^o®7oi^(n« + n(2)) (6.1) 
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with 



n(i) 
n(2) 




d^-^p 1 

{27T)D-ip\k\\q-k-p\' 
d^-^p [{q-k-p)kf 



(6.2) 



(6.3) 



(27r)^-iy (271)^-^21^131^-- A; 
By using dimensional regularization, Feynman-parametrization and usual integration for 
mulas we arrive at 

r(-Ke)r2(i-e) 



n(i) 



r2(i)(47r)i-T(2-2e)^ (2 



(6.4) 



and 

n(2) 



(4vr)2-(g2)-2- r2(l - 6)r(-l + 2e)T{\ - e) 
^ P(i)(47r)3 r(|-36) 

3r(f-e)r2(l-e)r(2e) 



4(4vr)3-2^(| - 2e)(-l + 2e)r(|)r(i)r(| - 3e) 
3r(l-e)r(2-e)(g2)i-2. 



+ 



2(47r)3-2T(|)r(|) 

r(26)r(|-6)r(2-6) 

(-l + 2e)r(f-e)r(f-3e) 



X 



(6.5) 



+ 



r(26) 

3-4e 



r(i-6) 

r(|-36) r(i-3e) ' r(|-36) 



4r(2 - e) ^ 4r(3 - e) 



2r(2e)r(| - e)r(2 - e) ^ 1(2 - e)r(l + 2e) 



r(|-e)r(|-3e) r(|-36) J- 

From eqs. (|6.4|) and (|6.5| ) the finite, renormalized contribution to the perturbation kernel 
may be extracted as 

ii4.a = 7o ®7o^^ln^. (6.6) 



Eq. ( |6.6D differs from eq. (|4.32|) by a simple factor proportional to a and therefore results 
in the mass shift 



8 l?TiQ!^ 

AM4 = ^ o (^i(^ + / + 1) + 7 + In ■ 



(6.7) 



The contribution from fig. 4.b is similar. After performing the integrations over the zero 
components we have 



n. 



ig^SN'Sab f d'^-^k f d'^-^p 



X 



-.2 



(27r)^-i 
2 



1 



p' 



+ 



(27r)^-i ^q-ky\p\\k\{q-p)^ 
{qp){pk){qk) 



(6.8) 
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This expression can be written in terms of integrals already solved in the course of the 
one loop calculation, and the outcome is 

Ho = -^^^|^^[I?^ + 2P(^-21n2)-2Plng2 (6.9) 

14 

+2 (7 - In 47r + 4 In 2) In g ^ + 2 In + const + 0(e)]. 

3 

Eq. (|6.9|) contains overlapping divergences and two graphs like 3. a with one or the other 
vertex replaced by a counterterm have to be added. After that, only an additive infinity 
survives which is subtracted by our usual renormalization prescription. Graph 4.b has a 
net contribution which is proportional to In^: 

n;?--^ln^|.. (6.10) 

The last two loop correction we are considering is the graph in fig. 4.c, whose contribution 
to the Coulomb gluon propagator can be written as 
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n. 



4.C 




rj 



{q - kYq 



1 



k^ 



+ 



(6.11) 



Eq. ( |6.11|) can be evaluated entirely in terms of gamma functions by a somewhat lengthy 
calculation, but following the same steps as above. The analytic result is 



n 



4.C 



2 Svr 



2^2 ' 



r2)i-2.(4^)2.r(e)r(2e)r2(2 - e)r(l - e)V{\ - 2e) 



r(|-26)r(|-3e)r(l + e) 



1 -2e)(l -4e). 

(6.12) 



Expanding in terms of e and properly renormalizing the result we finally obtain 



4.C 



'--InC + ln^'^ 



(6.13) 



487r4 ^ 3 /i- 

In table 3 we collect the results for the self energy parts of fig. 4. a to 4.c, apart from a 
factor -i^^q"^ 



graph 


n 


4.a 
4.b 


4 fl-' 

In^p 
flnC-iln'C 

3 /I^ 2 /i^ 


4.C 



Tab. 3 

The full calculation of the gluon self energy to two loops seems to be very involved in 
the Coulomb gauge. This is the more so, if massive fermions are included. However, in 
view of the results from the one loop calculation with massive fermions we may expect 
that nonvanishing masses tend to decrease the importance of such terms in practice to 
something that would be de facto numerically equivalent to O(a^). 

We try to circumvent these problems, for the time being, by the following argument, 
which also includes the three 'massless' quarks u,d,s. Because of the Ward identity for 
the Coulomb-vertex, it is clear from the theory of renormalization group that the same 
corrections can be obtained by expanding the running coupling constant with a two loop 
(gluons+u,d,s) input for the latter which provides also 'nonleading' logarithmic contribu- 
tions. Our present calculation in any case illustrates on the one hand the procedure to 
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be followed in a systematic BS perturbation theory. On the other hand, we believe that 
especially the computation methods for the notoriously difficult Coulomb gauge may be 
useful elsewhere. 

The beta function to two loops is renormalization scheme independent for massless 
quarks p8| and its two loop part has been calculated some time ago [p^] : 

(i{9) = -Pog' - Pig' - - (6.14) 

/^o = (6-15) 
1 38 

A = (4^^(102 --n,). (6.16) 
Here is the number of effective (massless) flavours and P{g) is the solution of 



In = / -— , 6.17 

-Jg Pig) 



which reads up to two loops 

1 



In 



/i 2po 



1 ^Pi^jHPo + Pig') 



(6.18) 



_g' g' Po g'iPo + Pig') 

Considering this as an equation for g = g{q'^) we 'undo' the renormalization group im- 
provement by expanding with 'small' g"^ oc a (cf. eq. (|2.6| ) ): 

33 -2n. g2 

m - 

167r fi 



a{q'^) = a <j 1 - In ^+ (6.19) 



Clearly the one loop term agrees with the calculation in subsect. 4.1 in the limit rrij —>■ 0. 
Nevertheless, this term should not be considered within the present argument because its 
error has the same order of magnitude as contributions of 0(a;^) . For the computation 

of the rest we need the expectation value of ^2 ■ This integral can be done analytically 
(Appendix A) and the result is: 

In^ — 2 

''^) = ;^{^;7 + ^2(n + / + l) + s„, + [^i(n + / + l)+7 + ln^]2}(6.20) 



with 



(f2 27™^ 12 am 



2(n-/-l)!"^2 (2/ + 1 + A;)! 



^ + 0' k\{n-l-l-ky 
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With eq. ( |6.2CI| ) we obtain for the mass shift, induced by the leading logs of the two loop 
vacuum polarization of the Coulomb gluon a contribution: 

^M^ioop = 2 o 127'[— + ^2(ri + / + l) + gnf + (^i(n + / + l)+7 + ln— )'] + 

1287r"^n^ [ 12 am 

+ 288(^i(r2 + / + l)+7 + ln^)l. (6.21) 

am J 

In this expression we have set the number of effective flavours equal to 3 as dictated by 
the number of sufficiently light quarks. Whether eq. ( p.21| ) really represents the full two 
loop quark- gluon vacuum polarization, numerically consistent with other terms 0(q;^), 
must still be checked in a calculation of the Coulomb gluon's self-energy to two loop order 
in the Coulomb gauge, i.e. going beyond the sample calculation here. 

6.2 QCD 2-loop Box Graphs 

It would be incorrect to extrapolate from the QED case the absence of corrections of 
0{a^), because gluon splitting allows new types of graphs. Our first example of a QCD 
box graph is fig. 5. a. Between the nonrelativistic projectors of the wave functions the 
perturbation kernel from this graph can be written effectively as 



-zH = -I2i 6 f J^lh (fco+p° + m)(g°-t°+pO 



m] 



(27r)4 (27r)4 [(^^ _ A;)2 - m2][(p2 - - m^]k \q- kft^iq- t)^ 
.(f-*H. + K£zMi^*WMzM|. (6.22) 



(t - k)' \ ' (t- kf 

After performing the integrations over the zero components and the external mo- 
menta can be set to (m, 0). This is justified a posteriori by the finiteness of the remaining 
terms. The resulting expression will thus only depend on q and m. The leading contribu- 
tion seems to be 

a^ 

H^.a oc — f^, (6.23) 

but a really reliable estimate or even a calculation of the coefficient of the leading term is 
not available yet. Supplementing the usual three powers of a from the wave functions for 
the computation of energy levels, we see that the graph 5. a indeed gives a contribution 
0(q;^). The qualitative result ( |6.23|) had been noted already in |jT6[, and |^. It should be 



noted, however, that also e.g. the graph 5.b may yield a contribution of the same structure. 
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A similar graph with crossed Coulomb lines (fig. 5.c) does not contribute because of the 
fact that the group theoretic factor vanishes. 

Now we investigate the 'X' graph in fig. 5.d for possible new contributions. As in the 
calculation of fig. 5. a it can be simplified to give 

Q f dH d'^k (A;o+p? + m)(t° + p°-m) 
—iH^d = oiq I - — —- — — — X 

The integration over k yields a l/\q\ divergence if q and p tend to zero. Contributions 
within the order of interest can therefore only come from possible poles after the t- 
integration. For simplicity we consider the part of eq. ( |6.24|) from the factor one in the 
second line: 

dH f +pl + m 



(27r)4[(t+p2)2-m2]t2(g_t)2 

-i /"i"^ yq^ — xp2 + p^ + m 



rl rl—x 

dx dy 
Jo Jo 



{Any Jo Jo {yq — xpY — yq^ — x{p^ — m 



2\ 



im . /"i-^ , X 



dx / dy — — + 0(a) 



{AtxY Jo Jo x'^m? + y{l — X — y)q 

In M + 0(a). (6.25) 



(47r)2m m 

Thus the part of graph 5.d, specified above, has a leading contribution of 

H - In 

32(47rj^|g]m m 

as g — >■ and therefore contributes to 0(a^ In a). The second part of graph 5.c gives a 
similar contribution, only with a different numerical factor. We conclude that - in contrast 
to the QED case [^]- box graphs may contribute to O(a^). As illustrated by the explicit 
calculations above, to O(a^) beside abelian QED type corrections a host of further 

non-abelian contributions can be foreseen. 

6.3 QED Correction 

As a rule, one is not forced to consider also electromagnetic effects in QCD calculations, 
but at high energies the strong coupling decreases, and in the case of toponium we expect 
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to be of the same order as aqEo- We can obtain this contribution by solving the BR 
equation for the sum of an QED and an QCD Coulomb gluon which results in the energy 
levels 



where Q is the electric charge of the heavy quark, that means 2/3 for toponium. Clearly 
even the 'leading' third term can only be separated from the effect of the second one to 
the extent that and aqEoif^) can be studied separately with sufficient precision. 

6.4 Weak Corrections 

While also weak interactions can usually be neglected in QCD calculations, this is not true 



in the high energy region. This is because the weak coupling scales like ^Gpvn? and this 
becomes comparable to the strong coupling if the fermion mass m is large. Even bound 



states through Higgs exchange may be conceivable Therefore we have to consider 
weak corrections and especially the exchange of a single Higgs or Z particle, assuming for 
simplicity the standard model with minimal Higgs sector. 
The Higgs boson gives rise to the kernel 

- lErnggs = -i\PiGprr? „ ^ ^ ^ i^Gpvr?—^^, (6.27) 

— mjj g 2 + 171^ 

where the notation should be obvious. 

Since we do not know the ratio am/niH which would allow some approximations if 
that ratio is small, we calculate explicitly the level shifts by transforming into coordinate 
space. As in Appendix A, we express the Laguerre polynomials in terms of differentiations 
of the generating function, do the integration and perform the differentiation afterwards 
to obtain an expression valid for arbitrary levels and Higgs boson masses: 

. , ^ Gp'm^Oi ^ , am , 

^MHiggs = -m ; J„K 6.28 

4V27r nmn 

with 



n2(l + a„)2" 



k=0 
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As an illustration some explicit results for the lowest levels are given in tab. 4. 



n 


1 


2 


2 


3 


3 


3 


I 








1 





1 


2 


Inlio-n) 


al 


al{2+al) 




a§(3+6a§+a|) 


a|(4+a§) 


4 


(l+ai)2 


4(1+02)4 


4(1+02)4 


9(1+^3)'' 


9(1+03)^ 


9(1+03)" 



Tab. 4 

It is evident that Eq. ( |6.28| ) will give a contribution of order Gprn'^a^ if the Higgs mass 
is comparable to the mass of the heavy quark and should therefore be taken into account 
in a consistent treatment of heavy quarkonia spectra to numerical order 0{a^). 

Next we consider the contributions of the neutral current, the single Z-exchange and 
Z- annihilation. 

For the Z-exchange we obtain 

Hz = -4V2GFml[Y{g, + g2i,)U> 'I ^ [Y{gi + 92!^%' p (6.30) 

q — rrii 

with 

9x = ^Ti -Qj sin^e^, (6.31) 

92 = -Iri, (6.32) 

where T3 is the eigenvalue of the diagonal SU(2) generator for the fermion /. If / is the 
top quark then = 1. 

Because mz is expected to be larger than am, it is possible to write the leading term 
of the mass shift resulting from (|6.30 ) in the form 



AMz = -4V2G^m|( ^'^^ ^^^'\ = (6-33) 

— 171^ 



m^^^{9l{?>-2S')-9l)5io 



V2' 



Tin 

where S is total spin of the quark-antiquark system. Therefore this expression gives rise 
to a singlet triplet splitting within the order of interest. 

In contrast to the gluon annihilation there is no colour structure involved in the Z- 
annihilation graph. Therfore this process does contribute. The corresponding energy shift 
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is easily evaluated: 



AMs=o = 7= :^Oio (6.34) 

V27r ri'^ 

2 2 

AMs=i = , ^Ms=o (6.35) 



Therefore this contribution also yields a singlet triplet sphtting. 



6.5 Schwinger-Christ-Lee Terms 



As mentioned in section nonlocal interactions have to be added to the Lagrangian 
of QCD in Coulomb gauge. We are not aware of any previous attempt to look whether 
these terms give contributions to bound state problems. In fact, by explicit calculation 
we found though that the Schwinger-Christ-Lee terms Vi and V2 in ( |2.1| ) only give rise 
to corrections to the propagator of the transverse gluon, and therefore can actually be 
neglected to 0(a'^). 



By analogy to the second ref. |T8[ we calculate the Vi term for SU(N) to 0{g^) 



J £r£r'£r"A^{f')Kij{f- f')Kjk{f- r")Al{f") (6.36) 



KM 



47r|p| 



%5(P) 



47r|p|' 



This corrects the gluon propagator by 

1 52 ■gAj^2 



Z[0]6J-{x,)54{x2) 16 



d'x I Mr'dh''--K,,K,kirMJ\- 



6 



5Ji 



(6.37) 



In momentum space 5G can be calculated by using dimensional regularization to give 

,4 Ar2 



SGtniQ.Q') = (2vr)S5'^''5(g-g')^^^(-5m„ + ^), 



(6.38) 



85 q2q2' 

which means that we have a mass shift with the same structure as the one transverse 
gluon exchange (cf. sect. (^), only suppressed by two more orders in a. 

Since the second term f 2 also represents a correction to the propagator of the transverse 
gluon, it can be estimated by the same method to contribute only in higher orders of a 
as well. 
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A Calculation of expectation values 

In sect. 4 and 6 we needed the expectation values of logarithmic potentials between 
Schrodinger wave functions. They can be obtained by 

if the expectation value (r^~^) is known analytically. This expression can be calculated 
by using the representation 

Lf+l ,ip) = hm , ^^^^ , (1 - zY^'^^e'^^ (A.2) 

of the Laguerre polynomials. This allows an easy evaluation of the integrations and the 
remaining differentiations can be done with some care: 



2n^-^ {n + iy. ' ' \n-l-l-k 



(A.3) 



Using now the Fourier transformations [pl|] 

7 + In /ir 



2 



In^ 



F[^] = (A.4) 

q ^ Inr 

In' 4 1 TT^ 

= — [- + 2(7 + lnH^] (A.5) 



one arrives immediately at eq.( 4.33 ) and (|6.20| ), respectively. 
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Figure Captions 
Fig.l: BS-equation for bound states 

Fig.2: Tree graphs (broken lines represent Coulomb gluons, curly lines depict transverse 
gluons, wavy lines represent a general gluon and solid lines stand for fermions). 

Fig. 3: One loop graphs and vertex corrections 

Fig.4: Two loop vacuum polarization 

Fig.5: Two loop box graphs 
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4.b 
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